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ABSTRACT 

In order to investigate the characteristics and influence of the magnetic field in evolved 
stars, we performed a follow-up investigatio n of our previous submillimeter analysis of the 
proto-planetary nebula (PPN) OH 231.8H-4.2 (ISabin et al.ll2014h . this time at 1.3 mm with the 
CARMA facility in polarisation mode for the purpose of a multi-scale analysis. OH 231 .8 h- 4.2 
was observed at ~ 2.5" resolution and we detected polarised emission above the 3cr thresh¬ 
old (with a mean polarisation fraction of 3.5%). The polarisation map indicates an overall 
organised magnetic field within the nebula. The main finding in this paper is the presence 
of a structure mostly compatible with an ordered toroidal component that is aligned with the 
PPN’s dark lane. We also present some alternative magnetic field configuration to explain the 
structure observed. These data complete our previous SMA submillimeter data for a better 
investigation and understanding of the magnetic field structure in OH 231.8H-4.2. 

Key words: Magnetic fields - polarisation - stars: AGB and post-AGB, stars individual: OH 
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1 INTRODUCTION 


The search for magnetic fields (hereafter B-fields) in evolved low 
and intermediate mass stars, from Asymptotic Giant Branch (AGB) 
stars to Planetary Nebulae (PNe), has been boosted in recent years 
by observations obtained thro ugh multiple techniques. T he detec¬ 
tions of s ynchrotron emission JPerez-Sanchez et al.|2013l) , Faraday 
rotati on ( Ransom et alj|2010h. maser emissi on ( Leal-Ferreira et al] 
I 2 OI 3 L IWol^ et al.l 2 OI 2 I: Arniri et ^ 1201 ih . dust continuum po- 
larisation 1 Sabin eLalJ |2007|), and mole cular line polarisation 
dGirart et all l2012l : IVlemmings et al.ll2012h have allowed us to de¬ 
termine the configuration and/or strength of the magnetic field Q. 
All these discoveries are leading to a better understanding of the 
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^ It is worth me ntioning that in P Ne, contrary to the case of AGB stars 
(see the work bv iLebre et al.|[20l3) . there is still no conspicuous and defi- 
nite measurements indicating the presence of a B-field at the ste llar surface 
jjordan et al j2012l ; lLeone et alj201 ill2014t[steffen et al.l2014ll . 


role magnetic fields play in the dynamics of evolved stars—in par¬ 
ticular, the role they play in shaping the circumstellar envelopes in 
PNe. 

In the case of dust grain polarisation, the assumption made is 
that in the presence of a magnetic field, non-spherical spinning dust 
grains will be aligned with their long axis perpendic ular to the B- 
field lines l lLazarianll2003l : Ikazarian & HoaindboOTl) . The inferred 
magnetic field orientations are therefore obtained by rotating the 
polarisation vectors by 90° allowing us to trace the direction of the 
magnetic field projected onto the plane of the sky. The polarisa¬ 
tion measurements at large wavelengths (i.e. from the far-infrared 
regime and above) are strongly correlated to the mass or intensity of 
the regions observed in the line of sight and thus the denser/heavier 
zones will have a larger weight on the determination of the polari¬ 
sation characteristics. _ 

Recently, our team dSabin et alj|20l4 hereafter Paper I) per¬ 
formed a new investigation of the dust continuum polarisation in 
the protoplanetary nebulae (PPNe) CRL 618 and OH 231.8+4.2 at 
345 GHz with the Submillimeter Array (SMA), an interferometer 
of eight 6m dishes. The array was used in compact configuration 
(~2" resolution). 
































2 L. Sabin et al. 


The main result of Paper I was the discovery in OH 231.8+4.2 
of an “X shaped” B-field pattern consistent with a poloidal configu¬ 
ration. There were also hints of a toroidal B-field component in the 
SMA data. The alignment between the magnetic field structure and 
the ^^CO(J=3 ^ 2) molecular outflow indicated a dynamically 
important field and the possible presence of a magnetic launch- 
ing mechanism c ompa tible with the magneto -centrifugal models of 
[Blackman et al.l l l200lh and iBlackmanI 1 20141) . 

In an effort to constrain these results we performed a follow¬ 
up investigation, this time at 1.3 mm, in order to obtain a multi¬ 
scale and multi-wavelength coverage of this source. The new po- 
larimetric observations provide more complete and spatially ex¬ 
tended magnetic field maps that will allow better comparisons with 
magneto-hydrodynamical (MHD) models. 

In §2 we describe our observations and data reduction, in §3 
we present the dust continuum and polarisation results for OH 
231.8+4.2. The discussion, which includes the comparison with the 
SMA data, and our conclusions are presented in §4 and §5, respec¬ 
tively. Finally, although this paper focused mainly on the source 
OH 231.8+4.2, we also chose to briefly present in the appendix sec¬ 
tion §A the non-conclusive CARMA detection for the other PPN 
studied in Paper I, namely CRL 618. 


2 OBSERVATIONS AND DATA REDUCTION 

The observations presented in this article were carried out with 
the Combined Array for Research in Millimeter-wave Astron¬ 
omy fCARMA. lBocdlrOO^ . The 15-element array (six 10.4-meter 
and nine 6.1-meter antennas) was used in D-configuration, which 
had baseline lengths between 11m and 148 m, and a resolution 
of ~2.5" at 230 GHz. The 1.3 mm dual-polarisation receivers 
jHull et alj|20llh were tuned to include both 1.3 mm continuum, 
the ^^CO(J = 2 —>■ 1) spectral line at 230.538 GHz, and the 
SiO( J = 5—^4) line at 217.105 GHz. We used a correlator config¬ 
uration with six 500 MHz wide bands to measure the dust contin¬ 
uum, and two 500 MHz wide band to map spectral-line emission. 
Our observations were performed in full-Stokes mode, which al¬ 
lowed all four polarisation crossproducts {LL, RR, LR, RL) to be 
measured in order to derive maps of Stokes 7, Q, U, and V. 

OH 231.8+4.2 was observed on 26 and 27 December 2013. 
The observing parameters for both sources are summarised in Table 
[T] It is worth noting the weather conditions during these track which 
was excellent for OH 231.8+4.2. 

The gain, passband, and flux calibrations, as well as the full 
processing up to the map production, were performed w ith the soft¬ 
ware MIRIAD dWright & Saulllll993l : ISault et alj20llh . The reduc¬ 
tion of CARMA polarisation data required two additional steps: 
XYphase and leakage calibrations. The XYphase calibration cor¬ 
rects for the phase difference between the L- and R-circular re¬ 
ceivers and is performed by observing a linearly polarised noise 
source with known position angle. The leakage calibration corrects 
any instrumental polarisation errors and is don e by observing a 
strong s ource over a range of parallactic angles. IHuU et aU 1 I2OI3L 
p.2) and [Hull etalJ l l2014l p.3) derived in their survey TADPOL 
(and for the same configuration used here) a typical band-averaged 
leakage of around 6%. More detailed information on the leakages 
can be found in these two references. The full description of polar- 
is ation observations with CARMA has been thoroughly discussed 
bv lHull et al] <20141) . 



Figure 1. HST Ha image of OH 231.8+4.2 (in logarithmic scale), with 
the superimposed CARMA 1.3 mm dust continuum emission (red contours 
with values of (0.010, 0.022, 0.040, 0.065, 0.100, 0.140, 0.180, 0.220, 
0.240) Jy/beam). 


Table 1. Summary of the CARMA observations 


Source 

OH 231.8+4.2 


Phase center 

0=07:42:17.0 

5=-14:42:50.2 

Obs. Date 

2013-12-27 

2013-12-26 

Gain calibrator 

0730-116 

0730-116 

Passband calibrator 

3C84 

3C84 

Flux calibrator 

Mars 

Mars 

Total project length (hrs) 

4.0 

4.0 

Time on source (hrs) 

2.37 

2.41 

Total Opacity! 

0.19 

0.16 


f Opacity at 230 GHz due to phase noise and atmospheric absorption 


3 POLARISATION AND MAGNETIC FIELD IN OH 
231.8+4.2 

3.1 Continuum polarisation 

The final continuum map of OH 231.8+4.2 has a synthesised beam 
size of ~4.6 x 2.4" with a position angle of 18°. The rms noise 
derived for the Stokes 7, Q, and U maps was aj = 1.7 mjy/beam 
and = 0.4 mJy/beam. The 1.3 mm continuum is elongated 
in the NE-SW direction (partially owing to the extended beam) 
and expands over ~ 15.3 x 8.8 arcsec^ (see Figure[T]l. We detected 
a peak intensity of ~240 mJy/beam centred on the coordinates 
a = 07‘'42“16!950, S = -14°42'49"95 and a mean flux of -70 
mJy/beam. 

Contrary to the “patchy” structure described in Paper I, 
the polarised continuum of OH 231.8+4.2 shown in Figure 
appears more uniform and extended (— 14.39x7.86 arcsec^). A 
comparison of both datasets is discussed in §4.1. We measured 
a peak polarised intensity (ppk) of about 3.6 mJy/beam and a 
mean of —1.7 mJy/beam (4.25 ct) over the whole polarised region. 
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Figure 2. Left Panel: Millimeter polarisation map of OH 231.8+4.2. The black contours showing the total dust emission are drawn in steps of 0.0241 Jy/beam 
X (0.2, 0.3, 0.5, 1, 2, 3, 4, 5, 6, 7, 8, 9). The grey scale image indicates the polarised intensity in Jy/beam. Finally, the polarisation vectors are drawn as orange 
line segments; the scale (bottom left) is set to 10%. Right Panel: Inferred magnetic field map obtained by rotating the original polarisation vectors by 90°. The 
map indicates a generally well organised magnetic field configuration within the whole structure. 



where we found the peak Ppk. The percentage polarisation (p) in 
this area reaches 12%, with a mean of 3.4%. Figure[2 shows that 
in this area the nine electric vectors (E-vectors) are well organised 
with a mean PA of —21 ± 7°. 

• The second area, or “lower distribution,” in the West is more ex¬ 
tended than the upper distribution but is similar in terms of p (peak 
at 12.3% and a mean of 3.5%). The main difference lies in the E- 
vector distribution. Figure [3 shows that while the region 2 has a 
mean PA of 4-34 ± 7° described by four E-vectors, region 3 has 
a mean PA of 4-14 ± 4° described by seven E-vectors, indicating 
different overall directions. The magnetic field appears to be well 
organised within each region. We point out that in both regions the 
rather large ppk is related to the outer edge of the nebula. 

We emphasise however, that the demarcation made between 
regions 2 and 3 is likely affected by the weight/number of E-vectors 
present in each area. Therefore instead of a clear distinction in the 
field orientation, we could actually be seeing a smooth variation of 
the magnetic field across the whole southern lobe. 


Figure 3. Zoomed-in view showing the variation of the polarisation position 
angle orientation. Three regions can be highlighted: region 1 with a mean 
PA of — 214i 7°, region 2 with a mean PA of 4-34 4i 7° and finally region 3 
with a mean PA of 4-14 ± 4°. Note that we need to add 90° to those num¬ 
bers to have the magnetic orientation angle. Combined with Figl^j-Right, we 
observe that while the first two zones are a reminiscent of the upper portion 
of the X-shape magnetic structure reported in Paper I, the third zone is a 
new and interesting element which fits with the “dark lane” or torus of OH 
231.8+4.2 (see §5). 


Based on the careful analysis of the polarised intensity, percentage 
polarisation, and position angle maps, the polarisation distribution 
could be separated in three regions within two main areas (Figure 

m-- 

• First an “upper distribution” in the East (region 1 in Figure O 


The drop of the fractional polarisation from the edges to the 
central area (near the total peak intensity), the so called “polar¬ 
isation hole,” is also seen in our map. The causes inducing this 
depolarisation effect have yet to be formally established, but are 
generally associated to three main phenomena. The first is beam 
depolarisation, where unresolved structures in the plane of the sky 
are averaged across the observing beam. Second, the low or null 
percentage polarisation can be linked to the poor alignment of the 
grains (i.e. low grain alignment efficiency) in high density regions. 
Finally, the fluctuation of the magnetic field along the line of sight 
is also likely to induce a change in the percentage polarisation. 

3.2 Magnetic field vs. ^^CO(J = 2 —> 1) molecular outfiow 

The establishment of a correlation between the distribution of 
magnetic fields and outflow orientations would help determine 


Jy/beam 
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Figure 4. Top: ^^CO(J = 2 —>• 1) molecular outflow detected with 
CARMA overlaid on the HST Ha map. The blue shifted lobe is shown as 
blue contours in steps of 8.137 Jy/beam ■ km/s x (-5, -4, -3, - 2 , -1, 1, 2, 3, 
4, 5, 6, 7, 8, 9); the red shifted lobe is shown as red contours in step of 5.586 
Jy/beam ■ km/s x (-5, -4, —3, -2, -1, 1, 2, 3, 4, 5, 6, 7, 8, 9). Bottom: Map 
of the magnetic field “vectors” superimposed on the ^^CO(J = 2 —>■ 1) 
molecular outflow. The grey scale image traces the dust continuum emis¬ 
sion. The best observed pattern is the globally perpendicular alignment of 
the southern vectors with the red-shifted outflow. 


whether PPN outflows are magnetically launched. Indeed, the 
co nditions for such launc hing mechanisms have been explained 
by [Blackman et alj (l2Q0lh . and rely on the presence of a poloidal 
field at small distance from the central engine and a toroidal field 
further away (see Paper I). However, we are still hindered by a lack 
of information on the strength of the magnetic field. 

The most direct path to establish this correlation 



Figure 5. Distribution of the thermal dust continuum emission at ~345 
GHz in blue (inner distribution-SMA data) and <^230 GHz in red (outer 
distribution-CARMA data), overlaid on optical H« data (HST). The combi¬ 
nation of high resolution data will allow us to perform accurate multi-scale 
magnetic analysis of the PPN (see §4). 

Bfield/outflow is to study the line polarisation (and trace the 
magne tic field) assuming the presence of a Goldreich-Kylafis 
effect dColdreich & KvlaM ll98lL Il982t) which predicts that the 
emission from rotating molecules can be polarised in the presence 
of a magnetic field. Th e strong ^^CO(J = 2—^1) line in 
OH 231.8-1-4.2 (see also lAlcolea et aljboOlh is therefore a good 
candidate to test the relationship between the magnetic field and 
the outflows. Unfortunately, the line polarisation mode is not 
available yet on CARMA. Therefore, as we did in Paper I, we rely 
on the comparison between the magnetic field derived from the 
continuum polarised emission and the direction and distribution of 
the CO emission line. 

Figure|4](top panel) shows the superposition of the CO emis¬ 
sion on the Ha HST image and we can immediately see the good 
correspondence between the blue-shifted emission and the north¬ 
ern optical outflow while the red-shifted CO outflow, roughly sim¬ 
ilar in size to its blue counterpart, only traces the interior (inner 
side) of the larger southern ionised outflow. In terms of polarisa¬ 
tion, the same Figure |4] (bottom panel), tends to indicate that the 
orientation of the “magnetic vectors” in the region 3 from Figure 
[3] is globally roughly perpendicular to the red-shifted (lower) CO 
outflow. The low number of “magnetic vectors” corresponding to 
the blue-shifted (upper) CO lobe prevents us from performing any 
more precise correlation between the magnetic field and the outflow 
direction. 

4 DISCUSSION 

The new CARMA millimeter data help to complete our understand¬ 
ing of the magnetic field configuration in OH 231.8-1-4.2. 
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4.1 CARMA/SMA comparative study 

The 1.3 mm data led to more extended maps of continuum intensity 
and polarisation. This allows us to perform a two-scale comparison 
of the magnetic structure (Figure [3. We find that the global mean 
continuum percentage polarisation of 3.4% is slightly lower than 
the 4.6% found with the SMA at -^345 GHz. 

In terms of configuration, the most interesting finding of this 
study is the detection of a group of organised “magnetic field 
vectors” globally perpendicular to both the ionised optical and 
molecular outflows; these orientations would then be coincident 
with an organised toroidal magnetic field. The presence of this 
structure is important for understanding how the magnetic field 
helps channel the outflow and how it affects the dynamics of OH 
231.8-1-4.2. 

By comparing the information of both the CARMA (with 
a maximum baseline of ~148m) and SMA (with a maximum 
baseline of ~77m) datasets and the respective polarisation maps, 
we observe that a toroidal magnetic structure, roughly aligned 
with the PPN equatorial plane, is (only) seen at larger scales. In 
Paper 1 we also suggested the presence of a toroidal magnetic 
structure, but the small set of vectors on which we relied (in the 
south-west area in the SMA data) were unlikely to belong to such 
configuration as the resulting torus would have an inclination in 
disagreement with the equatorial plane. Those vectors seen with 
the SMA are more likely to be part of the dipole field where the 
magnetic field lines are bending/curving. We also detected with 
CARMA what seems to be at first sight a separated poloidal 
component (north-east side) which could also be associated to a 
similar poloidal structure observed with the SMA in the northern 
lobes. However, the fact that those vectors are “only” seen on the 
eastern side combined with the possible variation of the vectors 
distribution in the South, stated in §3.1, might also suggest the 
presence of a single configuration in the form of a helical magnetic 
field, instead of two distinct field structure^ Proper modelling 
will be necessary to unveil the correct magnetic configuration. 

The detected magnetic configurations (and the scales on which 
they were detected) would be a first step in probing the presence of 
a magnetic launching me chanism in OH 231.8-1-4.2, as described 
bv iBlackman et alj l l200lh regarding proto-planetary nebulae. Fur¬ 
ther work, mostly modelling, is needed to fully confirm this phe¬ 
nomenon. 


4.2 Combined views of the magnetic field in OH 231.8-1-4.2 

OH 231.8-1-4.2 has been the subject of many p olarisation studies. 
We can cite first the work bv lEtoka et alj ( l2009l) . who investigated 
the configuration of the 1667 MHz OH maser emission and showed 
the presence a magnetic field aligned with the outflo w at a radius 
of~2” from the central star (matching the SMA scale. lSabin et al.l 
|2014|) . 

Then the hi gh resolution Very Long B aseline Array (VLBA) 
observations by iLeal-Ferreira et all ( 12012h . with a synthesised 
beam size of ~ 1.7 x ~0.9 mas, allowed the detection of 30 H 2 O 
masers. Only three spots exhibited linear polarisation, ranging from 
p=0.28% to p=1.15%. Circular polarisation was clearly detected 
in only one of the brightest maser spots, leading to a value of 


^ We note that in any case the helical configuration would be the combina¬ 
tion of both a poloidal and toroidal configurations. 


H| I = 44± 7 mG. In this case no configuration of the field could be 
accurately determined. The authors derived a surface stellar mag¬ 
netic field of -^2.5 G assuming a toroidal configuration. But if we 
assume that while going closer to the central star the field evolved 
from a mostly toroidal to a poloidal or dipole configuration (as our 
new results seems to suggest, but still keeping in mind the helical 
field hypothesis) this would then bring the stellar magnetic field 
strength to 140 G up to more than 8 kG respectively (L. Ferreira, 
private communication). A deeper magnetic map, would help con¬ 
firm this result. 


5 CONCLUSION 

We present a study of the dust polarisation in OH 231.8-1-4.2 at 1.3 
mm. The investigation follows up on a previous study conducted 
by our group at submillimeter range, and allows us to obtain a 
multi-scale map of the magnetic field distribution. In the case of 
OH 231.8-1-4.2, we successfully observed a globally organised 
magnetic field pattern in the form of a north-east poloidal section 
and what appears to be a toroidal configuration well distributed 
in the south. The presence of two distinct field configurations in 
the CARMA data is hampered by the geometry of the so-called 
poloidal field which is only coincident with one side of the outflow. 
An alternative explanation is the presence of a single configuration 
with a helical distribution including both the single-sided poloidal 
and the toroidal structures. 

The components found with CARMA complete those iden¬ 
tified with the SMA at smaller scales; thus, we conclude that the 
magnetic structure of OH 231.8-1-4.2 would consist of an inner 
poloidal field that is surrounded by an outer toroidal or possibly 
helical field . Further MHD modelling will be needed to properly 
describe the magnetic field configuration in OH 231.8-1-4.2 and link 
it (or not) to the signature of a magnetic collimation and launching 
mechanism. 

Much smaller scales (down to 0.5”) can now be reached with 
ALMA with its newly available polarisation system. The high 
resolution ALMA data will allow us to complete our multi-scale 
investigation of the magnetic field closer to the central engine, and 
to establish the role of the field in the late stages of stellar evolution. 
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APPENDIX A: THE PARTICULAR CASE OF CRL 618 

While this paper focused on the source OH 231.8+4.2, the proto¬ 
planetary nebula CRL 618 was also observed on 28 and 29 Decem¬ 
ber 2013. The observing parameters are summarised in Table lAl] 

When a source is very weakly polarised {Pf I ;$ 0.5%), very 
small changes to the leakage terms can cause drastic changes to the 
resulting polarisation position angles across the source. When we 
analysed the data for CRL 618, we saw position angles that var¬ 
ied by up to 90° from night to night depending on how exactly we 
reduced the data. The reason for this was that the source was ex¬ 
tremely bright at 1.3 mm (~2Jy/beam), but was very weakly po¬ 
larised (~ 0.5% polarisation at the intensity peak), which caused us 
to hit a dynamic range limit where the polarisation detected at very 
low levels can be caused by imperfections in the leakage solutions. 

The scatter in the real and imaginary parts of the leakage so¬ 
lutions on the two nights was ~0.01, which is larger than normal 
considering the 2 consecutive nights period; however, in this case 
those slight differences led to significantly different maps for Q and 
U. These small changes in leakages had caused virtually no vari- 
ations in the CARMA polarisation maps of protostellar cores (e.g. 
[Stephens et aT]|2013LlHull et al.ll20lll20Ti) , which were on aver¬ 
age at least a few percent polarised, and which tended to be much 
fainter sources where detection of polarisation was limited by the 
system temperature instead of by dynamic range. 

We therefore urge caution when interpreting CARMA obser¬ 
vations of sources with polarisation fractions of < 0.5%. While 
the very low polarisation percentage limit we find with CARMA is 
consistent with the value of 0.7% found with the SMA data, we are 
not able to confidently determine the B-field configuration. 


Table Al. CARMA observations of CRL 618 


Source 

CRL 618 


Phase center 

0=04:42:53.6 

5=36:06:53.4 

Obs. Date 

2013-12-28 

2013-12-29 

Gain calibrator 

3C111 

3C111 

Passband calibrator 

3C84 

0510+180 

Flux calibrator 

3C84 

Uranus 

Total project length (hrs) 

3.8 

5.5 

Time on source (hrs) 

2.43 

3.28 

Total Opacity! 

0.55 

0.86 


f Opacity at 230 GHz due to phase noise and atmospheric absorption 

This paper has been typeset from a TgX/ DTjX file prepared by the 
author. 









